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PREFACE 


The  investigation  reported  herein  was  autorized  by  the  Long  Beach 
Port  Authority  under  a contract.  Agreement  No.  WES  76-L,  dated  12  Sep- 
tember 1975. 

This  study  was  conducted  during  the  period  September  1975  to 
April  1976  in  the  Hydraulics  Laboratory  of  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  under  the  direction  of  Mr.  H.  B. 
Simmons,  Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  Chief  of 
the  Wave  Dynamics  Division  (WDD),  Mr.  D.  D.  Davidson,  Chief  of  the  Wave 
Research  Branch,  and  Mr.  C.  E.  Chatham,  Chief  of  the  Harbor  Wave  Action 
Branch.  Mr.  J.  R.  Houston  conducted  the  study  and  prepared  this  report. 

Messrs.  R.  R.  Bottin,  C.  Bryant,  and  A.  W.  Garcia  (WDD)  prepared 
the  finite  element  numerical  grids.  Mr.  D.  G.  Outlaw  (WDD)  developed 
the  computer  plotting  programs  and  aided  in  eliminating  errors  in  the 
finite  element  grids;  Mr.  Bottin  plotted  the  frequency  response  curves. 

Mr.  H.  L.  Butler  (WDD)  rewrote  the  matrix  solving  subroutine  of 
the  hybrid  finite  element  computer  program,  thus  reducing  significantly 
the  running  time  of  the  program. 

Drs.  H.  S.  Chen  and  C.  C.  Mei  of  the  Massachusetts  Institute  of 
Technology  provided  documentation  of  the  hybrid  finite  element  computer 
program  they  developed  and  materials  to  aid  in  its  utilization. 

A significant  proportion  of  the  numerical  computations  were  per- 
formed on  a CDC-7600  computer  at  the  Los  Alamos  Scientific  Laboratory 
(LASL),  Los  Alamos,  New  Mexico.  Mr.  P.  Johnson  and  Mrs.  Margaret 
McCormick  of  LASL  assisted  significantly  in  coordinating  visits  to  Los 
Alamos  and  scheduling  computer  time  on  the  LASL  system. 

Directors  of  WES  during  the  conduct  of  the  study  and  the  prepara- 
tion and  publication  of  this  report  were  COL  G.  H.  Hilt,  CE,  and 
COL  John  L.  Cannon,  CE.  Techincal  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary,  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)  units  as  follows: 


Multiply 


feet 

feet  per  second 

feet  per  second  per  second 

square  feet  per  second 


B Z 

0.3048 

0.3048 

0.3048 

0.09290304 


To  Obtain 

metres 

metres  per  second 

metres  per  second  per  second 

square  metres  per  second 


LONG  BEACH  HARBOR  NUMERICAL  ANALYSIS  OF  HARBOR  OSCILLATIONS 


EXISTING  CONDITIONS  AND  PROPOSED  IMPROVEMENTS 
PART  I:  INTRODUCTION 

Harbor  Oscillation  Problem 

1.  A harbor  excited  by  long  waves  with  a period  close  to  the 
natural  oscillation  period  of  the  harbor  experiences  a resonant  reaction 
that  may  produce  wave  heights  in  the  harbor  far  greater  than  the  in- 
cident wave  heights.  Generally,  the  vertical  movement  of  water  in 

the  harbor  is  still  only  a few  tenths  of  a foot  due  to  the  very  small 
heights  of  the  incident  long  waves.  However,  the  horizontal  movement 
of  water  may  cover  a distance  greater  than  10  ft*  as  a result  of  the 
long  wavelength  of  the  incident  waves.  Such  horizontal  water  movements 
cause  surge  and  sway  motion  of  moored  ships  that  may  hamper  unloading 
operations  and  damage  mooring  facilities. 

2.  A number  of  harbors  on  the  west  coast  of  the  United  States 
have  histories  of  ship  motion  problems  resulting  from  harbor  resonant 
oscillations.  It  is  likely  that  the  long-period  waves  responsible 
for  exciting  harbor  oscillations  permeate  the  Pacific  Basin.  Lack  of 
modern  harbor  development  may  prevent  the  occurrence  of  any  problem  for 
most  coastal  areas  of  the  Pacific.  Many  places  that  have  a problem 
may  merely  tolerate  it  through  ignorance  of  possible  remedial  measures 
with  the  result  that  the  problem  goes  unreported. 

3.  The  generating  mechanism  of  the  very  small  amplitude  long 
waves  that  excite  harbors  is  still  unknown.  However,  Munk's  surf-beat 
theory1  suggested  an  attractive  explanation  of  the  origin  of  these 
waves.  Surf  beats  result  from  the  superposition  of  two  trains  of  waves 
of  slightly  different  period.  The  beats  are  an  amplitude  modulation  of 

* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  is  presented  on  page  3. 
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a shorter  period  wave  train.  Swell,  for  example,  almost  always  exhibits 

an  irregular  beating  up  to  several  minutes  in  period. 

2 

U.  Carr  argued  that  local  surf  beat  could  not  cause  the  long- 
period  waves  known  to  excite  harbors.  Although  the  surf-beat  process 
produces  a long-period  vertical  water  motion,  the  really  important 
factor  of  long-period  horizontal  water  motion  is  lacking  in  surf  beats. 
Severe  harbor  surging  also  is  known  to  occur  during  otherwise  calm  seas 
with  identical  long-period  wave  activity  reported  along  hundreds  of 
miles  of  a coastline.  Carr  maintained  that  such  activity  is  not  com- 
patible with  long-period  wave  generation  by  local  surf  beat. 

5-  True  long-period  wave  trains  may  be  generated  from  surf  beats 
through  nonlinear  processes  in  the  surf  zone.  Munk^  estimated  that 
beaches  reflected  1 percent  of  the  energy  incident  upon  them  in  the  form 
of  these  long-period  waves.  It  is  perhaps  these  waves  reflected  from 
distant  shores  of  the  Pacific  which  cause  oscillation  problems  in  harbors 
on  the  west  coast  of  the  United  States.  Long-period  wave  activity  would 
then  increase  after  distant  storms  of  the  south  or  west  Pacific  Ocean 
created  large  swell  conditions  and  the  resulting  surf  beat  generated 

long-period  wave  trains  in  the  surf  zone  which  reflected  off  shorelines. 

3 

Wilson  has  noted  that  harbor  oscillation  in  Los  Angeles  and  Long  Beach 
Harbors  has  a seasonal  fluctuation  of  frequency  and  intensity  which  is 
greatest  during  the  summer  months  and  least  during  the  winter.  He  noted 
that  such  a pattern  is  almost  identical  with  that  found  for  Table  Bay 
Harbour,  Cape  Town,  South  Africa,  and  suggests  that  the  storms  in  the 
Southern  Hemisphere  are  responsible  for  the  long-period  waves. 

6.  Further  evidence  that  long -period  wave  trains  exist  from  time 

to  time  as  a general  ocean  condition  and  are  not  generated  by  local 

l* 

surf  beats  comes  from  a harbor  study  for  San  Nicholas  Bay,  Peru. 

It  was  concluded  from  horizontal  water  velocity  measurements  that 
long-period  waves  affecting  San  Nicholas  Harbor  were  progressive  waves 
unassociated  with  the  concurrent  swell  waves. 

Purpose  of  Study 

7.  The  purpose  of  this  study  was  to  investigate  harbor 
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oscillations  excited  by  long  waves  with  periods  from  1 to  10  min  using 
a finite  element  numerical  model.  Existing  conditions  of  Los  Angeles 
and  Long  Beach  Harbors  and  conditions  of  a proposed  harbor  expansion 
were  considered  to  ensure  that  the  expansion  would  not  cause  undesirable 
oscillations  in  existing  basins  or  characteristic  oscillations  of  its 
own  that  were  undesirable. 

8.  The  amplification  peaks  predicted  by  the  numerical  model  may 
be  much  larger  than  the  peaks  which  actually  occur  in  nature,  since  the 
model  neglects  all  dissipative  processes  except  radiation  of  energy  by 
a harbor.  However,  the  model  adequately  predicts  the  relative  severity 
of  various  modes  of  oscillation,  and  the  results  can  be  used  to  compare 
oscillation  characteristics  of  the  existing  Los  Angeles  and  Long  Beach 
Harbors  complex  and  proposed  harbor  expansions.  Information  about 
periods  of  resonant  peaks  and  the  spatial  patterns  of  wave  heights  and 
currents  when  areas  of  the  harbor  are  excited  also  will  be  used  to  de- 
termine both  the  incident  wave  period  ranges  which  need  to  be  tested 
and  the  optimum  gage  placement  during  future  wave  tests  in  the  Los 
Angeles  and  Long  Beach  Harbors  hydraulic  model. 


PART  II:  NUMERICAL  MODEL 

9.  Response  of  the  Los  Angeles  and  Long  Beach  Harbors  complex  to 
long  wave  excitation  was  determined  by  using  a hybrid  finite  element 
numerical  model  developed  recently  by  Chen  and  Mei  at  Massachusetts  In- 
stitute of  Technology.5  The  model  solves  the  following  generalized 
Helmholtz  equation: 

2 

V • [h(x,y) V(j>(x,y) ] + <i>(x,y)  = 0 (l) 

o 

where  <j>(x,y)*  is  the  velocity  potential  defined  by  U(x,y)  = 
with  U(x,y)  being  a two-dimensional  velocity  vector  and  w an  angu- 
lar frequency,  h(x,y)  is  the  water  depth,  and  g is  the  acceleration 
due  to  gravity.  Equation  1 governs  small  amplitude  undamped  oscilla- 
tions of  water  in  a basin  of  arbitrary  shape  and  variable  depth  forced 
by  periodic  long  waves.  It  has  been  further  assumed  that  the  flow  is 
irrotational . 

10.  The  Helmholtz  equation: 

2 

V (j>(x,y ) + ~ 4>(x,y)  = 0 ^ 

is  the  governing  equation  for  a constant-depth  ocean  region  outside  the 
basin. 

11.  For  a harbor  in  a semi-infinite  ocean  with  a straight  coast- 
line there  is  an  incident,  reflected,  and  scattered  wave.  The  scattered 
wave  has  a velocity  potential  <j>s  given  by 

00 

d>  = y a H (kr)  cos  n0  (3) 

vs  **  n n 
n =0 

where  a are  unknown  coefficients  and  H (kr ) are  Hankel  functions  of 
n 11 

the  first  kind  of  order  n . 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 

defined  in  the  Notation:  Appendix  A. 
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satisfies  the  radiation  condition  that  the  scattered 


wave  must  behave  as  an  outgoing  wave  at  infinity.  This  condition  is 
known  as  the  Sommerfeld  radiation  condition  and  may  be  expressed 
mathematically  as  follows: 


13.  Chen  and  Mei  used  a calculus  of  variations  approach  and 
obtained  a Euler -Lagrange  formulation  of  the  boundary  value  problem. 
The  following  functional  with  the  property  that  it  is  stationary  with 
respect  to  arbitrary  first  variations  of  <J>(x,y)  was  constructed  by 
Chen  and  Mei: 


F(4>)  = ff  1/2  {h(V<|>) 


where 

A = the  region  inside  the  harbor 
§ = line  integral 

<j>  = far  field  velocity  potential 

K 

<f>j.  = velocity  potential  of  the  incident  wave 

n = unit  normal  vector  outward  from  region  A 
a. 

a = boundary  of  region  A 

<j>  = total  velocity  potential  evaluated  on  boundary  a 

a 

lU.  Proof  was  given  by  Chen  and  Mei  that  the  stationarity  of 
this  functional  is  equivalent  to  the  original  boundary  value  problem. 

15.  The  integral  equation  obtained  from  extremizing  the  func- 
tional is  solved  by  utilizing  the  finite  element  method.  This  method 
is  a technique  of  numerical  approximation  that  involves  dividing  a 
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domain  into  a number  of  nonoverlapping  subdomains  which  are  called 
elements . 

16.  The  solution  of  the  problem  is  approximated  within  each 
element  by  suitable  interpolation  functions  in  terms  of  a finite  number 
of  unknown  parameters.  These  unknown  parameters  are  the  values  of  the 
field  variable  <p(x,y)  at  a finite  number  of  points  which  are  called 
nodes.  The  relations  for  individual  elements  are  combined  into  a 
system  of  equations  for  all  unknown  parameters. 

IT.  In  the  region  outside  the  basin,  the  velocity  potentials 
are  solved  analytically  in  terms  of  unknown  coefficients.  The  region 
is  considered  a single  element  with  an  "interpolation  function"  given 
by  Equation  3.  The  infinite  series  is  terminated  at  some  finite  value 
such  that  the  addition  of  further  terms  does  not  significantly  influence 
the  calculated  values  of  <j>(x,y)  . The  resulting  equation  is  combined 
with  the  system  of  equations  for  unknown  parameters  at  nodal  points 
within  the  basin  and  this  complete  system  is  solved  using  Gaussian 
elimination  matrix  methods. 

18.  n(x,y)  is  related  to  <(>(x,y)  through  the  linearized  dynamic 
free  surface  boundary  condition 


n(x,y) 


i 34>(x,y) 
g at 


(6) 


19*  The  horizontal  velocity  components  have  the  following  form: 


u(x,y ) 


£.  3n(x,y) 
w 9x 


v(x,y) 


-g  3n(x,y) 
w 3y 


(7) 


20.  The  hybrid  finite  element  method  (so  named  by  Chen  and  Mei 
because  the  method  involves  the  combination  of  analytical  and  finite 
element  numerical  solutions)  is  a steady-state  solution  of  the  boundary 
value  problem.  The  response  of  a harbor  to  an  arbitrary  forcing  func- 
tion can  be  easily  determined  within  the  framework  of  a linearized 
theory.  For  example,  an  arbitrary  incident  wave  amplitude  at  the 
harbor  mouth  in  the  absence  of  the  harbor  can  be  Fourier  decomposed  as 
follows : 
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f 
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bQ(t)  = J b(w)e  ivrt  dw  (8) 

_oc 

If  n(x,y,w)  is  the  response  amplitude  at  any  point  (x,y)  inside  the 
harbor  due  to  an  incident  plane  wave  of  unit  amplitude  and  frequency 
w , then  the  response  of  the  harbor  to  the  arbitrary  incident  wave 
amplitude  bQ(t)  is  given  by 


C(x,y,t)  = R 


<x 

/ 


b(w)n(x,y,w)e-iwt  dw 


(9) 


where  the  operation  R [ ] takes  the  real  part  of  the  bracketed 

quantity. 

21.  Therefore,  as  soon  as  n(x,y,w)  is  known  for  all  w , the 
harbor  response  to  an  arbitrary  incident  wave  can  be  calculated. 

22.  Chen  and  Mei^  verified  their  finite  element  numerical  model 
by  comparing  the  model's  calculations  with  analytical  and  experimental 
results  for  simple  wave  problems.  Figure  1 is  a comparison  of  analyt- 
ical and  numerical  results  of  wave  diffraction  by  a vertical  circular 
cylinder.  Similar  excellent  agreement  was  obtained  by  comparing  numeri- 
cal calculations  with  an  analytical  solution  for  a semicircular  harbor 
with  two  openings.  Experimental  data  for  the  oscillation  of  a rectangu- 
lar harbor  also  were  compared  with  the  finite  element  numerical  calcu- 
lations. Calculated  periods  of  the  resonant  peaks  agreed  with  experi- 
mental data;  however,  calculated  amplitudes  of  the  resonant  peaks  were 
larger  since  the  numerical  model  did  not  consider  frictional  effects. 

23.  Several  modifications  were  made  to  the  numerical  computer 
program  of  Reference  5 before  calculations  were  made  for  Los  Angeles 
and  Long  Beach  Harbors.  Slight  modifications  were  made  to  allow  the 
model  to  handle  harbors  with  variable  depths.  Chen  and  Mei^  applied 
their  numerical  model  only  to  constant  depth  problems.  The  subroutine 
which  solves  the  large  system  of  algebraic  equations  arising  from 
finite  elements  was  modified  to  be  more  efficient.  This  modification 
resulted  in  a decrease  in  the  computational  time  of  the  numerical  model 
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Figure  1.  Comparison  of  numerical  and 
analytical  results  for  wave  scattering 
off  a vertical  circular  cylinder^ 
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by  a factor  of  between  two  and  three. 

2k.  Velocity  calculations  also  were  added  to  the  numerical 
model. ^ Velocities  are  calculated  using  the  expressions  in  Equation 
The  slope  of  n may  be  obtained  at  any  point  (x,y)  within  an  element 
through  the  matrix  equation 


9n 

3x 

3n 

3y. 


G(n}e 


(10) 


where  the  matrix  vector 


(nl 


is  defined  as 


(11) 
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and  the  coordinates  of  the  three  nodes  (i,j,k)  forming  the  elements  are 
defined  as  (x^.y^),  and  (x^y^* 

25.  G is  the  element  slope  matrix  defined  as 


’aN. 

3N , 

9N, 

1 

iL 

k 

3x 

3x 

3x 

aN. 

3N . 

3N, 

1 

J. 

k 

_3y 

3y 

3y 

and  matrix  N is  known  as  the  "interpolation"  function  and  is  defined 
as 


N = {N. , N . , N } 
1 J k 


where 


”1  * ((Vj  - xjyk>  * (yk  - * <xj  - 

N.,N,  = similar  expressions  obtained  by  the  cyclical  permutation 
^ of  i ,j , and  k 

e 

A = area  of  the  element  A 

26.  The  {r)}6  matrix  is  known  from  the  finite  element  computa- 
tions. Together  with  the  interpolation  function  N , the  horizontal 
velocity  at  any  point  (x,y)  within  an  element  can  be  determined.  Ve- 
locities presented  in  this  report  were  calculated  at  element  centers. 

27.  The  finite  element  numerical  model  has  several  important 

advantages  compared  with  other  numerical  harbor  oscillation  models  such 

7 

as  that  of  Lee  and  Raichlen.  First,  the  model  allows  a variable  depth. 
Furthermore,  the  computational  time  requirement  for  the  model  is  sub- 
stantially less  than  computational  times  required  by  other  numerical 

Q 

models.  It  is  the  only  present  numerical  model  for  harbor  oscillations 
which  can  feasibly  calculate  resonance  effects  in  a harbor  as  large  as 
the  Los  Angeles  and  Long  Beach  Harbors  complex. 


PART  III:  APPLICATION  OF  NUMERICAL  MODEL 


Boundary  Conditions 


28.  The  boundary  condition  used  in  the  numerical  model  along  the 
shoreline  of  the  Los  Angeles  and  Long  Beach  Harbors  complex  is  that  the 
normal  component  of  the  velocity  be  equal  to  zero.  The  same  boundary 
condition  was  used  for  the  detached  breakwater  surrounding  the  pro- 
posed oil  tanker  terminal  at  Pier  J in  Long  Beach  Harbor.  The  break- 
water core  extends  above  still  water  and  is  impermeable  to  incident  wave 
energy  and  is  therefore  represented  as  a solid  barrier  in  this  study. 

The  eastern  entrance  to  the  tanker  terminal  basin  is  600  ft  wide  and  the 
basin  is  dredged  to  a depth  of  62  ft  mllw. 

29-  The  San  Pedro  breakwater  (Plate  l)  does  not  have  a solid 
core.  It  is  composed  of  class  B rock  with  an  average  diameter  of  ap- 
proximately 1.8  ft  and  has  a crown  at  mllw  of  fitted  concrete  blocks. 
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Keulegan  gives  the  following  formula  for  calculating  wave  transmission 
through  rock  structures: 


a^  = incident  wave  amplitude 
a^  = transmitted  wave  amplitude 
h = water  depth 
L = width  of  breakwater 
\ = wavelength  of  incident  waves 
T = wave  period 
Using  these  parameters 


30.  This  figure  is  a lower  bound  for  the  transmission  through  the 
San  Pedro  breakwater  for  a 3-min  wave  since  a^  = 0.05  ft  is  probably 
the  upperbound  wave  amplitude  for  a 3-min  wave  at  the  San  Pedro  break- 
water. This  upper  bound  can  be  deduced  from  prototype  measurements 
taken  in  Los  Angeles  and  Long  Beach  Harbors  by  WES.^  Long-period  wave 
activity  cannot  be  discerned  in  wave  gage  measurements  at  the  San  Pedro 
breakwater  because  of  the  large  wave  heights  in  the  swell  regime.  How- 
ever, a gage  off  Reservation  Point  in  Los  Angeles  Harbor  was  suf- 
ficiently protected  from  swell  that  visual  inspections  could  be  made  of 
long-period  wave  activity.  A total  of  107  visual  inspections  of  4-hr 
analog  charts  taken  during  ship  motion  occurrences  from  a year's 
prototype  recording  reveal  only  a single  instance  of  a wave  height  as 
great  as  0.1  ft  for  waves  with  periods  of  1 to  10  min.  The  single  re- 
cording was  of  a 6-min  wave  and  probably  results  from  the  large  6-min 
resonant  oscillation  that  occurs  in  the  East  Channel  of  Los  Angeles 
Harbor . 

31.  A 3-min  wave  was  considered  in  the  above  calculations  since 
longer  period  waves  of  the  same  amplitude  would  more  readily  be  trans- 
mitted through  the  breakwater  and  the  finite  element  grids  for  1-  to 
3-min  waves  (paragraphs  3^-36)  do  not  extend  as  far  as  this  breakwater. 

32.  Since  the  San  Pedro  breakwater  is  very  pervious  to  3-  to 
10-min  incident  waves,  it  is  not  represented  in  the  finite  element 
grids.  Although  the  breakwater  looks  impervious  from  the  water  as  a re- 
sult of  the  fitted  concrete  blocks  extending  up  from  the  waterline,  Wil- 
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son  noted  that  this  breakwater  has  long  been  considered  an  ineffective 


filter  against  long  waves  because  of  its  k3  percent  voids;  this  is  in- 
deed the  case. 

33.  The  Middle  breakwater  (Plate  l)  of  Los  Angeles  and  Long 
Beach  Harbors  has  a solid  core  to  -26  ft  mllw  which  is  covered  with 
a layer  of  class  B rock  and  a layer  of  class  A rock  (approximately  3 ft 
in  diameter).  The  rock  covering  of  the  solid  core  is  at  least  as  per- 
meable to  long  waves  as  the  San  Pedro  breakwater,  since  it  consists  of 
rock  the  same  size  or  larger.  Therefore,  the  Middle  breakwater  is 
represented  by  a depth  change  with  the  core  impermeable,  and  small 
finite  element  triangles  are  placed  along  the  breakwater  in  the  numer- 
ical grid  to  properly  represent  it. 


Numerical  Grids 


3**.  The  numerical  grids  used  in  the  numerical  model  computations 
are  shown  in  Plates  3-12.  The  elements  are  all  triangular  and  often 
equilateral.  The  lengths  of  the  triangle  sides  were  in  general  selected 
so  that  they  were  one-eighth  of  the  local  wavelength  or  less.  Chen  and 
Mei'’  compared  the  finite  element  solution  of  diffraction  by  a vertical 
cylinder  with  an  exact  analytical  solution  and  found  a global  error  of 
0.7  percent  when  they  used  an  element  length  of  one-tenth  the  incident 
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wavelength.  They  defined  a global  error  as  — (local  error),  where 
the  summation  is  performed  over  n finite  element  nodes.  The  local 
error  is  defined  as 


where  n is  che  magnitude  calculated  by  the  finite  element  method  at 
a nodal  point  and  n is  the  theoretical  value.  In  the  present  study, 
it  was  judged  that  element  lengths  equal  to  one-eighth  of  the  local 
wavelength  or  less  would  provide  sufficient  accuracy.  Element  lengths 
vary  in  the  grids  in  relation  to  local  depths.  Element  lengths  near 
harbor  entrances  are  made  much  smaller  than  may  be  strictly  necessary 
in  order  to  ensure  proper  propagation  of  waves  into  inner  basins. 


35-  Five  finite  element  grids  were  used  for  each  case  tested  and 
are  tabulated  below.  It  was  not  feasible  to  use  a single  grid  to  cal- 
culate harbor  response  for  incident  waves  with  periods  from  1 to  10  min. 
For  example,  1-  to  3-min  grids  could  not  be  constructed  to  cover  the 
entire  harbor  complex  because  of  computer  memory  and  running  time  limi- 
tations. Separate  grids  were  used  for  3-  to  4-,  it-  to  6-,  and  6-  to 
10-min  calculations  to  reduce  the  cost  of  operating  the  model.  The 
computational  time  of  the  model  is  approximately  proportional  to  the 
number  of  nodes  in  the  grid  times  the  bandwidth  squared.  The  bandwidth 
is  the  greater  of  either  the  maximum  difference  for  all  elements  be- 
tween two  adjacent  node  numbers  in  an  element  of  the  grid  or  the  sum  of 
the  number  of  nodes  on  the  grid  semicircle  plus  the  number  of  terms 
retained  in  truncating  the  infinite  series  of  Equation  3.  In  general, 
the  bandwidth  increases  approximately  linearly  with  a decrease  in  the 
lengths  of  the  triangular  element  sides.  The  number  of  nodes  increases 
approximately  as  the  inverse  square  of  the  triangle  side  lengths.  Thus, 
a grid  for  3-min  waves  (eight  nodes  per  wavelength)  requires  approxi- 
mately l6  times  more  computational  time  than  one  for  6-min  waves  (eight 
nodes  per  wavelength) . 


Existing  Conditions  Revised  Conditions 


Grid 

1 

1 

to 

2 min 

Grid  1A 

1 

to 

2 min 

Grid 

2 

2 

to 

3 min 

Grid  2A 
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to 

3 min 

Grid 

3 
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to 

4 min 

Grid  3A 
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to 

4 min 
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k 
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to 

6 min 

Grid  4a 

4 

to 

6 min 

Grid 

5 

6 

to 

10  min 

Grid  5A 

6 

to 

10  min 

3 6.  Many  of  the  numerical  grids  for  existing  conditions  were 
designed  so  that  they  could  be  more  readily  modified  to  construct  grids 
for  the  harbor  expansion  conditions.  However,  the  design  of  the  harbor 
was  modified  several  times  during  the  course  of  this  study.  Therefore, 
some  element  patterns  in  the  grids  may  have  no  particular  purpose. 

Grids  also  were  constructed  to  minimize  the  bandwidth,  and  this  require- 
ment occasionally  governed  triangle  shapes  for  areas  of  the  grids. 
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PART  IV:  RESULTS 


Comparison  of  Prototype  Data  and  Numerical  Calculations 
for  Existing  Conditions 


37-  Prototype  wave  gage  recordings  were  collected  at  Los  Angeles 
and  Long  Beach  Harbors  from  26  May  1971  to  29  June  1972  by  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station.^  Energy  spectra  were 
determined  for  10  time  periods  of  the  year  during  which  ship  motion 
problems  were  reported.  Spectral  analyses  were  performed  via  Fast 
Fourier  Transform  (FFT)  techniques.  Variability  of  spectral  estimates 
was  decreased  by  smoothing  over  10  adjacent  frequency  bands.  Ten  bands 
were  chosen  so  that  there  would  still  be  sufficient  frequency  resolution 
to  define  the  long-period  energy  for  wave  periods  as  great  as  10  min. 

The  data  record  length  analyzed  for  each  spectral  analysis  was  24,576 
sec  (6.83  hr)  and  the  lag  cime  between  successive  record  start  times  was 
chosen  to  be  3 hr.  Further  details  of  the  analyses  can  be  found  in 
Reference  10. 

38.  Tables  1-5  give  estimates  of  wave  periods  and  spectral 
energy  densities  of  peak  energies  for  prototype  wave  gages  5,  6,  7*  8, 
and  l4;  Plate  1 shows  the  locations  of  these  gages.  The  starting  times 
in  Tables  1-5  give  the  year,  Julian  day,  and  hour  on  a 24-hr  clock. 

Each  of  the  10  time  periods  lasts  approximately  2 days.  Data  do  not 
exist  for  some  time  periods  as  a result  of  gage  downtime.  Gage  lU  was 
in  operation  only  during  one  time  period. 

39.  Tables  1-5  give  the  wave  period  with  the  greatest  spectral 
energy  density  within  a range  of  periods.  The  range  of  periods  is  such 
that  the  spectral  lines  within  the  range  cannot  be  resolved  by  the 
spectral  analysis;  for  example,  in  the  range  480-360  sec,  the  adjacent 
frequency  lines  corresponding  to  periods  of  459.36  and  387.02  sec  cannot 
be  resolved  by  the  spectral  analysis  and  should  be  considered  the  same. 

40.  Plate  13  is  a time-dependent  spectral  plot  of  gage  5 data 
for  the  time  period  of  1500  November  27-0300  November  29,  1971. 

The  variability  of  energy  levels  in  this  plate  illustrates  the 


nonstationarity  of  the  variance  in  energy  bands.  Long-period  waves  are 
generally  transient  (time-dependent)  phenomena,  and  some  caution  must  be 
observed  in  the  application  of  spectral  analysis  techniques  to  nonsta- 
tionary processes.  However,  Tables  1-5  do  give  an  indication  of  char- 
acteristic energy  bands  for  each  gage  location. 

41.  It  should  be  noted  that  the  finite  element  nodal  calculations 
provide  only  the  response  characteristics  of  Los  Angeles  and  Long  Beach 
Harbors.  This  response  must  be  multiplied  by  the  input  spectrum  of 
waves  for  a comparison  with  prototype  data.  However,  the  input  wave 
spectrum  is  not  known  since  long  waves  arriving  from  the  deep  ocean 

had  heights  too  small  to  be  measured  by  gages  near  the  outer  breakwaters 
during  the  prototype  collection.  Still,  a correspondence  exists  be- 
tween resonant  peaks  calculated  by  the  numerical  model  and  those  mea- 
sured in  the  prototype.  Plates  14-38  are  plots  of  the  wave-height 
amplification  factor  calculated  by  the  finite  element  numerical  model 
for  prototype  gages  5,  6,  7,  8,  and  l4.  The  frequency  response  for  the 
existing  conditions  of  Los  Angeles  and  Long  Beach  Harbors  shown  in 
Plate  1 is  compared  with  the  frequency  response  for  the  revised  condi- 
tions shown  in  Plate  2.  The  wave-height  amplification  factor  is  defined 
at  each  point  as  the  wave  height  at  the  point  divided  by  twice  the  in- 
cident wave  height.  This  definition  of  amplification  factor  is  tradi- 
tional and  is  a result  of  the  fact  rhat  the  standing  wave  height  for  a 
straight  coast  with  no  harbor  would  be  twice  the  incident  wave  height 
due  to  the  superposition  of  the  incident  and  reflected  waves. 

42.  Table  1 indicates  that  there  is  no  resonant  peak  for  gage  5 
data  between  8 and  10  min.  The  numerical  calculations  presented  in 
Plate  18  also  show  no  resonant  peak  for  this  period  range.  Table  1 
shows  that  large  peaks  occur  at  387  and  237  sec.  Plates  l6-l8  show 
large  peaks  at  396  sec  and  222  sec  and  no  peaks  from  240-360  sec. 

Neither  the  numerical  computations  in  Plate  15  nor  the  prototype  data  of 
Table  1 indicate  peaks  from  120-210  sec.  The  prototype  data  have  peaks 
at  100,  93,  and  64  sec.  The  numerical  model  calculates  peaks  at  92,  89, 
and  63  sec  as  seen  in  Plate  l4.  Excellent  agreement  on  the  resonant 
peaks  is  apparent  between  prototype  data  and  the  numerical  model 
calculations. 
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43.  Table  2 shows  a large  resonant  peak  for  gage  6 at  387  sec, 
small  peaks  at  237  and  l6o  sec,  fairly  large  peaks  at  93  and  86  sec, 
a moderate  peak  at  80  sec,  and  small  peaks  in  the  range  65-75  sec. 

Plates  19-23  show  a large  peak  at  396  sec,  small  peaks  at  222  sec  and 
158  sec,  a fairly  large  peak  at  88  sec,  moderate  peaks  at  80  and  92  sec, 
and  a small  peak  at  66  sec.  Therefore,  resonant  peaks  calculated  by  the 
numerical  model  and  measured  in  the  prototype  are  in  excellent  agree- 
ment for  gage  6. 

44.  Gage  7 has  large  resonant  peaks  in  Table  3 at  564  and  334 
sec,  fairly  large  peaks  at  2l6  and  150  sec,  and  a small  peak  at  JO  sec. 
Plates  24-28  show  the  frequency  response  calculated  by  the  numerical 
model  for  gage  J.  Plate  28  indicates  a moderate  peak  at  518  sec  and 
part  of  what  may  be  a large  peak  beyond  600  sec.  These  peaks  probably 
correspond  to  the  564-sec  peak  of  Table  3,  since  the  spectral  lines 
adjacent  to  the  564-sec  line  are  located  at  733  and  459  sec;  therefore, 
the  564-sec  peak  contains  energy  over  a large  range  of  periods.  The 
numerical  model  calculates  a large  broad  peak  from  300  to  350  sec  with 
a maximum  at  338  sec  (Plate  27);  a large  sharp  peak  at  280  sec  is  also 
shown  in  Plate  27.  Such  a sharp  peak  may  not  be  important  in  the  pro- 
totype since  there  would  be  little  incident  energy  concentrated  within 
such  a small  range  of  wave  periods.  Part  of  the  energy  in  this  resonant 
peak  also  may  appear  in  the  334-sec  spectral  line  in  the  prototype  data 
analysis.  The  numerical  model  also  predicts  a moderate  peak  at  245  sec, 
a small  peak  at  216  sec  (Plate  26),  and  a large  peak  at  177  sec  (Plate 
25)  which  probably  corresponds  to  the  150-sec  peak  in  the  spectral 
analysis.  Plate  24  shows  a fairly  large  but  sharp  peak  at  74  sec. 

This  peak  may  correspond  to  the  small  peak  at  70  sec  in  the  spectral 
analysis.  Again,  the  sharpness  of  the  peak  would  limit  the  occurrence 
of  the  oscillation  in  the  prototype  data. 

45.  Gage  8 displays  fairly  large  resonant  peaks  in  Table  4 at 
334  and  l60  sec,  moderate  peaks  at  564,  190,  120,  and  85  to  100  sec, 
and  a small  peak  at  76  sec.  Plate  33  shows  part  of  a peak  which  may 
occur  beyond  600  sec  and  correspond  to  the  564-sec  peak  in  the  spectral 
analysis.  Plate  32  shows  a fairly  large  broad  peak  at  338  sec  and  a 
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very  sharp  peak  at  280  sec , which  again  may  not  be  important  in  the 
prototype  as  a result  of  the  sharpness  of  the  peak.  There  is  a large 
peak  in  Plate  31  at  2l8  sec  which  corresponds  to  the  199-sec  peak  in  the 
prototype  data.  Plate  29  shows  peaks  at  ll6,  90  to  105,  78  to  90,  and 
60  to  75  sec.  These  peaks  are  similar  to  those  noted  in  Table  4.  The 
very  sharp  peak  at  108  sec  is  probably  not  important  in  the  prototype. 

46.  Table  5 gives  estimates  of  wave  periods  and  spectral  energy 
densities  of  peak  energies  for  one  time  period  for  gage  l4.  Wave 
energies  during  this  time  period  were  fairly  low  (Tables  1-4).  Resonant 
peaks  in  Table  5 occur  at  387,  237,  160,  100,  90,  70,  and  64  sec. 

Plates  34-38  show  peaks  at  396,  222,  l60,  100,  92,  and  from  65  to  80  sec. 
The  periods  of  resonant  peaks  calculated  by  the  numerical  model  and  mea- 
sured in  the  prototype  are  in  excellent  agreement  at  gage  l4. 

Comparison  of  Harbor  Response  for  Existing 
and  Revised  Conditions 


47.  Plates  14-38  show  a comparison  of  plots  of  the  wave-height 
amplification  factor  as  a function  of  wave  period  at  gages  5,  6,  7,  8, 
and  14  for  existing  and  revised  harbor  conditions.  Plates  1 and  2 show 
the  existing  harbor  and  the  proposed  modifications.  The  modifications 
include  landfills,  dredge  cuts,  and  a small  breakwater.  Wave  periods 
from  1 to  10  min  are  considered. 

48.  Plates  39-63  show  a comparison  of  plots  of  the  normalized 
maximum  current  velocity  versus  wave  period  at  gages  5,  6,  7,  8,  and 
l4  for  existing  and  revised  harbor  conditions.  The  normalized  maximum 
current  velocity  is  a maximum  velocity  over  one  wave  period,  and  it  is 
normalized  so  that  multiplication  by  the  incident  wave  amplitude  in 
feet  gives  velocity  in  units  of  feet  per  second. 

49.  Plates  64-76  show  a comparison  of  plots  of  the  wave-height 
amplification  factor  versus  period  for  three  sites  in  Long  Beach  Harbor, 
sta  LB-1,  LB-2,  and  LB-3  (Plate  l).  Prototype  wave  gages  are  not  lo- 
cated at  these  sites;  however,  the  three  sites  may  experience  resonance 
problems.  Plates  77-89  show  a comparison  of  plots  of  the  normalized 


maximum  current  velocity  versus  period  for  these  same  three  sites. 

50.  Plates  90-97  show  plots  of  the  wave-height  amplification 
versus  period  for  two  sites  (sta  LB-4,  LB-5)  off  Pier  J in  the  area 
where  the  oil  terminal  is  proposed  (Plate  l).  Plates  98-105  show  plots 
of  the  normalized  maximum  current  velocity  versus  period  for  the  same 
two  sites.  These  two  sites  were  chosen  to  illustrate  existing  condi- 
tions in  the  area  where  the  oil  terminal  is  to  be  located.  The  plots 
indicate  that  there  are  no  large  resonant  peaks  between  2 and  10  min 

at  these  sites. 

51.  Plates  106-114  show  plots  of  the  wave-height  amplification 
factor  versus  period  for  the  three  tanker  terminals  in  the  oil  terminal 
of  the  revised  Long  Beach  Harbor.  Plates  115-123  show  plots  of  the 
normalized  maximum  current  velocity  versus  period  for  the  same  three 
tanker  terminals.  The  main  resonant  peak  for  all  three  terminals 
occurs  at  a period  of  220  sec.  The  wave-height  amplification  factor 

is  similar  for  all  three  terminals  at  this  period.  This  is  an  indication 
that  at  this  period  the  basin  between  the  detached  breakwater  and  the 
Pier  J landfill  has  an  oscillation  mode  in  which  the  water  surface  in 
the  center  part  of  the  basin  rises  and  falls  fairly  uniformly.  This  is 
indeed  the  case  as  can  be  seen  in  Plate  146.  Plate  l46  shows  contours 
of  the  wave-height  amplification  factor  over  the  entire  grid  area.  The 
basin  is  experiencing  a "pumping  mode"  type  oscillation  in  which  the 
water-surface  level  inside  the  basin  rises  and  falls  almost  uniformly. 

52.  The  resonant  response  of  the  oil  terminal  basin  to  long 
wave  excitation  shown  in  Plates  106-123  is  very  similar  to  the  response 
calculated  in  a previous  report‘d  which  used  the  same  numerical  model 
but  a much  finer  numerical  grid  that  only  covered  the  basin  and  the 
immediate  surrounding  area.  The  numerical  model  predicts  an  amplifica- 
tion factor  peak  of  approximately  4 at  a 4-min  period  for  both  grids; 
the  resonant  response  is  smoother  (fewer  small  peaks)  for  the  finer 
grid,  however.  This  difference  probably  results  from  the  fact  that  the 
finer  grid  covers  only  a small  region  around  the  terminal  area.  There 
are  no  external  influences  for  this  case  such  as  nearby  resonances 
which  may  couplb  with  the  oil  terminal  basin  response  or  reflected 


waves  from  other  parts  of  the  harbor. 

53-  Plates  1^-89  indicate  that  there  is  no  radical  difference  be- 
tween the  resonant  oscillations  characteristic  of  the  present  and  re- 
vised harbors.  The  amplification  factor  and  normalized  velocity  plots 
have  similar  shapes  with  resonant  peaks  located  at  corresponding  periods 
The  resonant  peaks  for  the  revised  harbor  conditions  are  in  general 
shifted  to  slightly  lower  periods  and  have  small  amplitudes,  for  example 
Plates  16,  18,  and  23. 

5^.  There  are  some  very  sharp  resonant  peaks,  for  example,  the 
108-sec  resonant  peak  shown  in  Plate  29,  which  occur  in  response  plots 
for  one  harbor  configuration  but  not  another.  This  difference  is  prob- 
ably attributable  to  difficulties  in  locating  the  peaks  and  not  dif- 
ferences in  response.  Response  calculations  were  initially  made  every 
2,  3,  and  5 sec  for  incident  waves  with  periods  from  1 to  3,  3 to  7, 
and  7 to  10  min,  respectively.  Resonant  peaks  were  then  further  defined 
by  calculation  intervals  of  0.5  to  2 sec  and  sometimes  0.1  to  0.25  sec. 
If  a resonant  peak  appeared  for  one  harbor  configuration  and  not  another 
there  was  a search  for  the  missing  peak.  Occasionally  it  was  difficult 
to  locate  a very  sharp  peak.  These  peaks  are  probably  not  important 
in  the  prototype  since  they  are  so  sharp  that  only  small  amounts  of 
energy  in  the  incoming  wave  spectrum  would  be  amplified. 

55-  Plates  12h-l6l  are  plots  of  contours  of  wave-height  amplifi- 
cation and  plots  of  normalized  maximum  current  velocity  for  many  of  the 
resonant  peaks  shown  in  Plates  lU-89.  The  velocities  are  represented 
by  arrows  whose  centers  lie  at  element  centroids.  Water  particles  move 
horizontally  back  and  forth  in  a direction  parallel  to  the  arrows.  The 
amplitude  of  the  incident  waves  is  assumed  to  be  equal  to  1 ft.  The 
velocities  are  in  units  of  feet  per  second  for  this  amplitude  incident 
wave.  Incident  waves  in  the  prototype  may  have  amplitude  orders  of 
magnitude  smaller  than  1 ft  and  the  velocities  they  produce  are  cor- 
respondingly smaller. 

56.  Plates  12U-l6l  basically  illustrate  patterns  of  wave-height 
amplification  and  current  velocity  and  show  how  the  harbor  is  oscillat- 
ing. Quite  often  the  resonances  are  localized.  For  example,  the  U-ir.in 


oscillation  shown  in  Plate  144  is  localized  in  the  gage  5 channel  and 
does  not  extend  throughout  the  Southeast  Basin.  This  oscillation  is 
the  fundamental  mode  for  the  channel  with  an  antinode  forming  at  the 
end  of  the  channel  and  a node  at  the  mouth.  It  is  interesting  to  note 
that  a simple  Merian's  formula  calculation  also  predicts  a fundamental 
period  of  oscillation  of  the  channel  equal  to  222  sec.  The  depth  of 
the  channel  is  assumed  to  he  58  ft  and  the  length  2U00  ft.  This  simple 
calculation  works,  of  course,  because  this  channel  has  a very  simple 
shape. 

57-  Plate  156  is  an  illustration  of  a more  complicated  oscilla- 
tion which  involves  all  of  Southeast  Basin.  A node  cuts  the  basin  in 
half.  Antinodes  which  are  out  of  phase  develop  at  the  ends  of  basin  6 
and  the  gage  5 channel.  Plate  156  indicates  that  the  amplitude  of  the 
396.5-sec  oscillation  should  be  approximately  the  same  at  prototype 
gages  5 and  6.  The  square  root  of  the  ratio  of  energies  in  the  360- 
to  U80-sec  period  range  for  gages  5 and  6 (Tables  1 and  2)  is  very  close 
to  unity  for  the  10  time  periods  of  the  prototype  data.  The  wave  energy 
is  proportioned  to  the  square  of  the  wave  amplitude.  This  is  an  addi- 
tional example  of  the  close  agreement  between  prototype  data  and  calcu- 
lations from  the  numerical  model. 

58.  Plates  124-l6l  aid  in  understanding  the  resonant  oscilla- 
tions which  develop  in  the  Los  Angeles  and  Long  Beach  Harbors  complex. 
Knowledge  of  oscillation  patterns  would  be  helpful  in  determining  re- 
medial actions  to  limit  oscillations,  the  proper  location  of  wave  gages 
in  a hydraulic  model  of  the  harbors,  and  areas  of  low  surge  (antinodal 
locations) . 
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PART  V:  CONCLUSIONS 


59-  The  good  agreement  between  prototype  measurements  in  the  Los 
Angeles  and  Long  Beach  Harbors  complex  and  calculations  of  the  hybrid 
finite  element  numerical  model  confirms  the  reliability  of  this  numeri- 
cal model  as  a viable  engineering  tool  for  studying  resonant  oscilla- 
tions of  complex  harbors.  Amplitudes  of  resonant  peaks  predicted  by  the 
numerical  model  will  always  be  too  large  since  the  model  does  not  include 
any  dissipative  effects  other  than  radiation  dissipation.  However,  the 
model  can  predict  relative  severities  and  is  therefore  extremely  use- 
ful in  comparing  the  oscillation  characteristics  for  existing  condi- 
tions with  those  for  various  proposed  modifications.  The  rapid  calcu- 
lation time  of  the  numerical  model  makes  it  feasible  to  apply  the  model 
to  large  harbors  and  broad  incident  wave  period  bands.  This  numerical 
model  can  be  used  as  an  aid  to  hydraulic  model  tests  since  it  pinpoints 
the  periods  of  resonant  peaks  and  shows  oscillation  patterns  which  are 
useful  in  optimizing  wave  gage  placement  to  measure  extreme  conditions. 

60.  The  numerical  model  calculations  indicate  that  resonant 
oscillations  in  existing  harbor  facilities  will  not  be  radically  altered 
by  the  proposed  modifications  of  the  harbors  complex.  The  resonant 
oscillations  are,  in  general,  somewhat  smaller  for  the  modified  condi- 
tions of  the  harbor  complex  than  for  existing  conditions.  There  are 
some  peaks  which  are  larger  for  the  modified  conditions;  however,  they 
are  not  drastically  larger. 

61.  There  is  only  one  well-defined  resonant  peak  from  1 to  10  min 
for  the  proposed  oil  terminal  basin  modification.  The  basin  tested  has 
a 600- ft  eastern  breakwater  opening  and  a 62-ft  depth.  The  resonant 
peak  has  a period  of  approximately  U min  and  an  amplification  factor 

of  approximately  h.  This  amplification  factor  is  relatively  small 
compared  with  amplification  factors  for  many  resonant  peaks  shown  in 
Plates  lU— 89 • Whether  or  not  this  peak  could  cause  significant  ship 
motion  problems  depends  upon  the  natural  resonant  modes  of  the  moored 
tankers  berthing  in  the  basin.  There  are  no  resonant  peaks  in  the  oil 
terminal  basin  of  such  a magnitude  that  they  would  obviously  preclude 


safe  berthing  at  the  terminals  given  proper  mooring  practices  and 
adequate  mooring  equipment. 

62.  Detailed  information  on  moored  ship  response  as  a function 
of  incident  wave  amplitude  and  frequency,  the  incident  wave  spectrum, 
and  the  resonant  response  of  the  harbors  complex  to  incident  wave  ex- 
citement must  be  known  before  definitive  conclusions  can  be  made  on 
the  precise  degrees  of  ship  motion  in  Los  Angeles  and  Long  Beach  Har- 
bors. Numerical  model  calculations  presented  in  this  report  indicate 
that  the  resonant  response  of  the  harbors  complex  would  probably  not 
be  changed  enough  to  cause  an  increase  in  the  severity  of  resonant 
oscillations.  In  general,  the  resonant  oscillations  for  the  proposed 
modifications  are  somewhat  smaller  than  those  for  existing  conditions 
The  few  peaks  which  increased  in  magnitude  were  very  sharp  (narrow 
frequency  band)  and  probably  will  not  result  in  a significant  problem 
due  to  their  limited  bandwidth. 
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Table  1 


:!it imates  of  Wave  Periods  and  Spectral  Ene;  ry  Densities 
,-f  Peak  Knergies  for  Gage  5 


Ranges  oi  Wave  Periods 
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of  Time  Hist  or: 
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71-168-0300 


71-289-1600 


71-329-0300 


* Wave  periods  in  seconds. 

*•  Spectral  energy  densities  in  10‘ 


Table  2 


estimates  of  Wave  Periods  and  Spectral  Energy  Densities 
of  P»Jak  Energies  for  Gage  6 


;es  of  Wave  Periods 
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* Wave  periods  In  seconds. 

'*  Spectral  energy  densities  in  10" 


Starting  Time 

: 

of  Time  History 

- aMMI 

!i?0-3^0 

’ . 


Estimates  of  Wave  Periods  and  Spectral  Energy  Densities 


;es  of  Wave  Periods 


Starting  Time 


of  Time  Histo] 
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71-31 u-1500 


71-316-0900 


80-60 

69.52 

iu.23 

69.52 

iu.ea 

69.52 
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69.  52 

3.69U 

69.52 

18.82 

69.52 

7.0U9 
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71-289-1600 


71-304-1200 
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71-329-0300 


72-155-0600 


Wave  periods  in  seconds. 
Spectral  energy  densities  In  10' 
No  data  for  this  time  period. 


Estimates  of  Wave  Periods  and  Spectral  Energy  Densities 


of  Peak  Energies  for  Gage  8 


* Wave  periods  in  seconds. 

**  Spectral  energy  densities  in  10‘ 
t No  data  for  this  time  period. 


-2  2 
i ft  -sec. 
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Table  5 


an 


Estimates  of  Wave  Periods  and  Spectral  Energy 
Densities  of  Peak  Energies  for  Gage  lU 


Starting  Time 

Ranges 

of  Wave 

Periods 

of  Time  History 

1800-U80 

1+ 80- 360 

270-210 

170-150 

139-95 

72-155-0600 

733.61* 

387.02 

237-1+5 

160.01 

100.93 

186.9** 

13.59 

22.28 

91.33 

1+6.85 

90-80 

80-65 

65-60 

51-1+6 

72-155-0600 

89.86 

69.52 

6b.  08 

1+9.80 

11.80 

22.7  b 

18.21 

20.50 
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Boundary  of  region 

Amplitude  of  incident  wave 

Amplitude  of  transmitted  wave 

Region  inside  harbor 

Amplitude  of  frequency  component  w 

Incident  wave  amplitude 

Average  rock  diameter 

2 

Acceleration  due  to  gravity,  32.2  ft/sec 
Element  slope  matrix 
Water  depth,  ft 

Hankel  function  of  the  first  kind  of  order  n 
Imaginary  number 
Wave  number 
Integer 

Interpolation  function 

Unit  normal  vector  outward  from  region  A 

Spherical  coordinate,  ft 

Real  number 

Time,  sec 

Wave  period 

Velocity  in  x-direction,  ft/sec 
Total  horizontal  velocity,  ft/sec 
Velocity  in  y-direction,  ft/sec 
Angular  frequency,  radians/sec 
Cartesian  coordinate,  ft 
Cartesian  coordinate,  ft 
Unknown  coefficient 
Coefficient 
Area  of  element 
Gradient  operator,  ft-1 
Spherical  coordinate,  radians 
Wavelength  of  incident  waves 
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£ Response  of  harbor 

2 

<j>  Total  velocity  potential,  ft  /sec 

2 

(>  Total  velocity  potential  evaluated  on  boundary  a , ft  /sec 

a 2 
(>T  Velocity  potential  of  incident  wave,  ft  /sec 

1 2 

Far  field  velocity  potential,  ft  /sec 

” 2 
|>  Scattered  wave  velocity  potential,  ft  /sec 
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